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1. Introduction 

The present paper addresses the following problem. Let w{x) be a classical discrete 
weight function, and consider discrete symplectic and orthogonal ensembles associated 
to these weights, see Sections [2] and [3] for precise definitions. By methods of Random 
Matrix Theory (see e.g. Tracy and Widom [JTj) we find m-point correlations in terms of 
Pfafiians with 2x2 matrix kernels. The problem is to express these kernels in terms of 
the orthogonal polynomials associated with our weight function. 

It is well known that Pfaffian expressions including 2x2 matrix kernels appear in 
analysis of the orthogonal and the symplectic ensembles of Random Matrix Theory, see, 
for example, Mehta p6], Forrester [1^, Tracy and Widom [H], Borodin and Strahov [9]. 
Such kernels also play a role in the works on the crossover between matrix symmetries, see 
Pandey and Mehta [37], Mehta and Pandey [39], Nagao and Forrester [32], superpositions 
of matrix ensembles, see Forrester and Rains [211 [22], ci-nd also in the multi- matrix models 
combining matrices of different symmetries, see Nagao [271 [29] . Forrester and Nagao 
[18] , Borodin and Sinclair [8J give 2x2 matrix kernels for ensembles of asymmetric real 
matrices, in particular, for the eigenvalue statistics of Real Ginibre Ensemble. Vicious 
random walkers, random involutions, Pfaffian Schur process are examples of problems 
from combinatorics and statistical physics where Pfaffian formulas including 2x2 matrix 
kernels arise, see Nagao and Forrester [33], Nagao, Katori and Tanemura [M], Nagao [28] . 
Forrester, Nagao and Rains [19j, Borodin and Rains [4j, Vuletic [12]. 

Often these kernels are constructed in terms of skew-orthogonal polynomials. Then 
a question arises how to compute the skew-orthogonal polynomials, and how to handle 
the Christoffel-Darboux sums involving them. In some cases explicit formulas for skew- 
orthogonal polynomials have been given in terms of related orthogonal polynomials, and 
matrix kernels and their asymptotic values have been computed. This approach is devel- 
oped in Nagao and Wadati [53], Brezin and Neuberger [TD], Adler, Forrester, Nagao and 
P. van Moerbeke [T], see also Forrester [T7], Chapter 5. Nagao [2S] provides the list of 
the cases when the expressions of skew orthogonal polynomials in terms of the classical 
orthogonal polynomials are explicitly known, see Table 1. 

Our definitions of discrete symplectic and orthogonal ensembles in Sections |2] and [3] 
are motivated by relations with 2;-measures on Young diagrams with the Jack parameter 
^ = 2, as it is described in Section jH For ensembles obtained in Section H] skew-orthogonal 
polynomials are not known, and we use a discrete version of the method developed by 
Widom in ^3] in the context of orthogonal and symplectic ensembles of Hermitian matri- 
ces. Widom [12] gives general formulas expressing entries of 2 x 2 matrix kernels in terms 
of the scalar kernels for the corresponding unitary ensembles. Whenever the logarithmic 
derivative of the weight in the definition of orthogonal or symplectic ensemble under con- 
siderations is a rational function, the entries of the 2x2 matrix kernels are expressible 
in terms of orthogonal polynomials, and are equal to the scalar kernel plus extra terms. 
Similar results for ensembles with Laguerre-type weights were obtained in physical liter- 
ature, see Sener and Verbaarschot [38], Klein and Verbaarschot [24J. These papers show 
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that the number of extra terms is finite, which leads to universahty of correlation kernels 
for such ensembles. 

Formulae obtained in Widom [13] are especially convenient for the asymptotic analysis 
since the asymptotics of polynomials associated to rather general classes of weights is 
known, see Deift, Kriecherbauer, McLaughlin, Venakides, and Zhou [121 [13], Bleher and 
Its [2]. This enables one to use Widom' s formulae in the proof of the universality for 
the orthogonal and symplectic ensembles, see Deift and Gioev [HI [15], Deift, Gioev, 
Kriecherbauer, and Vanlessen [TB], Stojanovic [40j . 



i ne weignt 
defining the families 


The orthogonality 

set 


The family 

of the orthogonal 

polynomials 


w{x) = e~^^ 


(-00, +oo) 


Hermite polynomials 
(see [Il[ini[35]) 


w{x) = x°'e~^ 


(0,+cx)) 


Laguerre polynomials 
(see [Il[30]) 


w{x) = (1 -x)"(l - xf 


(-1,1) 


Jacobi polynomials 
(see P [301 [35]) 


w{x)=[{l+x)\{l-xy]-' 
L is an even integer 


Z 


Hahn polynomials 
(see [^) 


w{x) = 1 


{0,1,. ..,L} 


Hahn polynomials 
(see [2^) 


w{x) = 


Z>o 


Meixner polynomials 
M„(x;c= l,g) (see [19]) 



Table 1. The cases in which skew-orthogonal polynomials are explicitly known. 



Our results for discrete symplectic and orthogonal ensembles are of the similar kind as 
those obtained in Widom ^3]. Whenever a discrete analog of the logarithmic derivative 
of the weight is a rational function the matrix kernels are expressible in terms of the or- 
thogonal polynomials associated to weight in the definition of the ensemble. This is used 
to work out the cases of the Charlier ensemble (w(x) = ^, x G Z>o), and the Meixner 
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ensemble {w{x) = ^^^c^, x G Z>o). As an application, we compute the continuous limit 
of our formulas corresponding to the degeneration of the Meixner orthogonal polynomials 
to the Laguerre orthogonal polynomials. 

Acknowledgements. We thank Percy Deift for discussions. This work was supported 
in part by BSF grant BSF-2006333. The first named author (A.B.) is grateful to Dim- 
itri Gioev for helpful discussion during early stages of this work. He was also partially 
supported by the NSF grant DMS-0707163. 

2. Main results for discrete symplectic ensembles 

Let w{x) be a strictly positive real valued function defined on Z>o with finite moments, 
i.e. the series X]i;ez>o '^{x)^^ converges for all j = 0, 1, 

Definition 2.1. The A'^-point discrete symplectic ensemble with the weight function w 
and the phase space Z>o is the random A"-point configuration in Z>o such that the prob- 
ability of a particular configuration xi < . . . < xjv is given by 

N 

Prob {xi, . . . , xn} = Z]^\ ^^w{xi) J]^ {xi - Xjf{xi - Xj - l){xi - Xj + 1). 

i=l l<i<j<N 

Here Zjv4 is a normalization constant which is assumed to be finite. 

In what follows Z^^ is referred to as the partition function of the discrete symplectic 
ensemble under considerations. 

Introduce a collection {P„(C)}^o '^^ complex polynomials which is the collection of 
orthogonal polynomials associated to the weight function w, and to the orthogonality set 
Z>o. Thus 

• P„ is a polynomial of degree n for all n = 1, 2, . . ., and Pq = const. 

• If m 7^ n, then 

J2 Prn{x)Pn{xMx) = 0. 

xe'z>o 

For each n — 0,1,... set (pn{x) = {Pn, Pn)Z^^'^ Pn{x)w^^'^{x), where (.,.)«, denotes the 
following inner product on the space C[Q of all complex polynomials: 

(/(C), ^(0).:= E mgi^M^)- 

a;eZ>() 

We call i^n the normalized functions associated to the orthogonal polynomials P„. 
Let 7i be the space spanned by the functions (po, (pi, . . .. 
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Definition 2.2. Suppose that there is a 2 x 2 matrix valued kernel Kn4:{x, y), x,y & Z>0; 
such that for a general finitely supported function rj defined on Z>o we have 

N 

(1 + ?7(xi)) J]^ {xi- XjY{xi- Xj -l){xi- Xj + 1) 

(xi<...<XAr)CZ>o j=l l<i<j<N 

= Vdet {I + 7]Kn4), 

where Kna, is the operator associated to the kernel KN4{x,y), and r] is the operator of 
multiplication by the function rj. K^iix^y) is called the correlation kernel of the discrete 
symplectic ensemble defined by the weight function w{x) on the phase space Z>o. 

An explanation of the term "correlation kernel" can be found in Tracy and Widom |41j . 
§2, 3. 

We introduce the operators D^,D^ and e which act on the elements of the space 7i. 
The first and the second operators, and D^, are defined by the expression: 

{D^f){x)= J2 D±ix,y)fiy), 
where the kernels D± (x, y) are given explicitly by 



w[x] 

(2.1) D+{x,y) = y\— — — x,yeZ>o, 

wix + 1) 



wix — 1 ) 

2.2) D_{x,y) = ^1 ^^^^^ ' 6,,i,y, x,yeZ>o- 

The third operator, e, is defined by the formula 



+ 00 



V- / w(2m) wi2m + l)w{2m + 3)...w{2k + l) 

^ y w[2k + 1) w[2m)w[2m + 2) . . . w[2k) 

(2.3) 

... , ^ / w{2k) w{2k + l)w{2k + ?>) ...w{2m + l) 

(eo9) (2m + 1) = > \ — —— , , , i ) r -(p{2k), 

^^'^ ' ^Yw(2m+1) w{2k)w{2k + 2) ...w{2m) 

where m = 0, 1, . . .. 

To make sure that etp is well defined for any y9 G 7i, we impose an additional assumption 
on the weight function w: we assume that 

w{x — 1) di{x) 



x> 1 

w{x) d2{x)^ ~ 

for some polynomials d\ and such that deg d\ > deg d2 and if deg di = deg d2 then 
\im di{x) / d2{x) > 1. This implies, in particular, that w{x — l)/w{x) > const > 1 for 
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X ^ 1, and one easily verifies that the series defining {eip){x) converges for any ip E Ti 
and X G Z>o. 

Let us also introduce the operator S'7V4 which acts in the same space Ti, and whose 
kernel is SN4ix,y). To write down S^Aix^y) exphcitly, introduce 2N x 2N matrix M^^"* 
whose j, k entry (j, A; = 0, 1, ... , 2N — 1) is 



(2.4) 



2;fci2,>0 



where D := Dj^ — D^. 

Proposition 2.3. The matrix M^^^ is invertible. 
All the proofs are delayed until Section [51 

Write {M^'^^)~^ = (jJ^fj^), and define the kernel SN4{x,y) by the formula: 

2N-1 

(2.5) SN4{x,y) = <pj{x)nfi^<pk{y), 

j,k=0 

where x,y E Z>o. 

Theorem 2.4. The operator K^^ (see Definition \2.2\] is expressible as 



K 



s 



SnaD. 



Remark 2.5. As it is clear from the proof, the operator K^a for the discrete symplectic 
ensemble can be also represented as 

V+S'Ar4 — V+S'Ar4V_ 



(2.6) 



K 



NA 



In the formula just written above the operators V+, V_ are defined by 

(V+/) (x) = (/(x + 1) - f{x)) , 



(V-/) {x) 



wix + 1) 



'wix — \] 



w[x) 



(/(x)-/(x-l)) 



Let T^AT be the subspace of Ti spanned by the functions (/jq, V^i, 
Km the projection operator onto TIn- Its kernel is 

2N-1 

KN{x,y) = E (pkix)(pkiy). 



ip2N-i- Denote by 



fc=0 



It is convenient to enlarge the domains of D and e, and to consider the operators 
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e: n + DH^n + eH. 
It is not hard to check that these operators are mutual inverse. Denote by -D-^jv 
restriction of the operator D to Hn- 

Theorem 2.6. The following operator identity holds true 

Dhn^na = {Ihn+dHn ~ [-0' KN\KNe) ^ K^. 

The next Theorem gives the condition on the weight function w{x) under which the 
operator Sna can be written in an exphcit form. 

Theorem 2.7. Let w{x) be a weight function such that for x > 1 

w{x — 1) di{x) 
w{x) d2{xy 

where di,d2 are polynomials of degree at most m satisfying the assumption above, and 
di{0) = 0, ^2(0) ^ 0. Then 

n 

(2.7) [D,KM]KN = J2^^^^^, 

i=l 

where a®h denotes the operator with the kernel a{x)b{y), ipi, . . . ,ipn elements of 
Hn, and ipi, . . . ,4^ 

n are elements of T~ij^. Assume in addition that the matrix T^j — 
5ij + {eipi, ipj), i, j = 1, . . . ,n is invertible. Then 

n 

(2.8) SMi = eKM-J2 iT-%{ei,i) ® (K^e^^,) . 

Set d2{x) = const -(x — ai)"^ . . .{x — aiY^ , and let be the order of '^^('^^^ at oo. 

As will be clear from the proof of Theorem 12.71 the number n from (12.71) is bounded by 

I 

+ Yl ^ai- III Proposition 110. ll we show that n = 1 implies T = 1. 

i=l 

Corollary 2.8. // the commutation relation between the operators D and K^q takes the 
form 

[D, Kn] = \{ipi ® + ® V^i), 
where ipi G Ti.N,4'2 ^ H'^n, and A is some constant, then 

Sn4: = — Xeip2 ® e^pi. 

The general formalism described above can be applied in particular to the Meixner 
and to the Charlier symplectic ensembles. The weight function for the Meixner symplec- 
tic ensemble is by definition the weight function associated to the Meixner orthogonal 
polynomials: 

(2.9) WMeixner{x) = ■^^c'', X G Z>o, 
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where /? is a strictly positive real parameter, and < c < 1. The weight function of the 
Charlier symplectic ensemble is defined by 

(2.10) Wcharlier{x) = —r, X G Z>o, 

x\ 

where a > 0. wchariier is the weight function defining the classical Charlier orthogonal 
polynomial^. 

Theorem 2.9. a) If w{x) is the Meixner weight with the parameters c and (3 defined by 
equation l{2.9\) . then the operator Sna whose kernel is defined by equation l{2.5\) takes the 
following form: 

J2N(2N + 3-1) ^ , , , 
(1 -c)Vc 

where the operator Kn has the kernel 



.01 IN f \ ^2Nc{2N + /3 - 1) 'f2N{.x)^2N-l{.y) - V2N-l{.x)^2N{.y) 

(2.11) KN[x,y) = , 

1 — c X — y 

the functions {v^fc(a;)}^o ^'^^ normalized functions associated to the Meixner orthog- 
onal polynomials, the operator e acts by the formula 



^ V (^ + m)i (m+ 2)/+! 



(e^) {2m + l) = V-cf: J [ ^"i; ^(2- - 21), 

V i-rn-^)i+i 

where m = 0, 1, . . and the functions ipi, ip2 o^^e defined by the expressions 

(2.12) ^i(x) = v^^^ - v/2iV + /?-l^Hi^, 

X + p a; + p 

(2.13) Ux) = v/2iV + /5-l^^ - v/2A^^^^^. 

X + p — 1 X + p — 1 

b)If w{x) is the Charlier weight with the parameter a (see equation Ii2.10\) ). then the 
operator 8^4 whose kernel is defined by equation 112. 5\) takes the following form: 



2N 

Sn4 = (^Kn + \ {eip2N) ® {eip2N-l) , 

V a 

where the operator Kn has the kernel 

-(f2N{x)(f2N-l{y) - f2N-lix)(f2N{y) 



KN{x,y) = -VWg 



x-y 



^FoT definitions and basic properties of the classical discrete orthogonal polynomials sec Ismail 
and also Koekoek and Swarttouw [Ml 
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the functions {^k{x)}'^^Q are the normalized functions associated to the Charlier orthog- 
onal polynomials, and the operator e acts as follows: 

V 2 ^ (m+ 

(ey.) (2m + 1) = a/I E / ¥^(2"^ ' 20- 

V 2 ^ (-m- 

3. Main results for discrete orthogonal ensembles 

Definition 3.1. The 2A^-point discrete orthogonal ensemble with the weight function 
W and the phase space Z>o is the random 2A^-point configuration in Z>o such that the 
probability of a particular configuration xi < . . . < X2n is given by 

Prob{xi, . . .,X2n} = 

2N 

^Ni n ^i^i) n i^j ~ if ~ i^ odd for any and Xi is even, 

i=l l<i<i<2N 

0, otherwise. 

Here Z^i is a normalization constant. 

In what follows we assume that the weight function W{x) is such that 

(3.1) W{x - l)W{x) = w{x), for x > 1, and W{0) = w{0), 

where w{x) is a strictly positive real valued function on Z>o satisfying the same conditions 
as the weight function in the definition of the discrete symplectic ensemble in Section |2l 

Definition 3.2. Suppose that there is a 2 x 2 matrix valued kernel Kni{x, y), x,y E Z>o, 
such that for an arbitrary finitely supported function r] defined on Z>o we have 



2N 

Ww{x,){l + r]{xi)) JJ {xj-Xi) = Vdet {I + r]Km) 

{xi<...<X2NKZ>a «=1 1<*<J<2A' 



where is the operator associated with the kernel KNi{x,y), and rj is the operator of 
multiplication by the function rj. Then the kernel of Kj^i is called the correlation kernel 
of the discrete orthogonal ensemble defined by the weight function W{x) on the phase 
space Z>Q. 

As in the symplectic case, details on the correlation functions can be found in Tracy 
and Widom [tl], §2, 3. 

Let e and D be as in the previous section, and let w{x) in the definitions of these 
operators be given in terms of the weight function W{x) by formula (13.11) . Introduce 
the operator Sm which acts in the same space H, and whose kernel is S]\fi{x,y). To 
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write down S]\fi{x, y) explicitly introduce 2N x 2N matrix M'^^^ whose j, k entry (j, k 
0,1,...,2A^-1) is 

(3.2) Mf^= J2 <x,y)v,{x)My)- 

Proposition 3.3. The matrix M^^^ is invertible. 

Write (M^-^^)^^ = {jJ^fl), and define the kernel SNi{x,y) by the formula: 

2Af-l 

(3.3) SNi{x,y) = (^j(x)/iJVfc(y), 

j,k=0 

where x,y E Z>o. 

Theorem 3.4. The operator K^i (see Definition \3.2\] is expressible as 



K 



Nl 



Sni^ Sni 



Denote by e-j^^ the restriction of the operator e to Tijv. Recall that is the projection 
operator on T^at. 

Theorem 3.5. The following operator identity holds true 



Theorem 3.6. Let w(x) be as in Theorem\2. 71 Then 



n 



(3.4) [e,Kj,]Kj, = Y, 



i=l 



where rji,...,rjn are elements o/Hn, and fii,...,rin are elements ofHj^. Assume in 
addition that the matrix Uij = 6ij + {Dr]i,fjj), i, j = 1, . . . ,n, is invertible. Then 

n 

(3.5) Smi = DK^ - J2iU-%{Dm) ® [K^Dr],) . 

The number of terms in fl3.4p is bounded in the same way as n in Theorem 12. 7[ 

Corollary 3.7. // the commutation relation between the operators D and Kjsr takes the 
form 

[D, Kn] = Xi^Ji ® ^2 + ^"2 ® V'l), 

where ipi € T^Af? ^"2 ^ '^'nj and A is some constant, then 

Sni = DKn - )\ip2 ® i>i- 



DISCRETE ENSEMBLES 11 

Theorem 3.8. For Meixner or Charlier weight (see equations 1^2. 10\) and 1^2. 9\) ), we have 
the following expressions for the operator S^i : 

a) For the Meixner orthogonal ensemble 

(3.6) Sm = DKn H -r. ^—^ i)2 ® ^i, 

(1 -c)Vc 

where the functions ipi and ip2 are defined by equations ^2.12\) and 112. 13\) respectively. 

b) For the Charlier orthogonal ensemble 

(3.7) 5'Ari = DKn H V2N ® '^2N-i- 

a 

4. Discrete symplectic and orthogonal ensembles related with 

^-measures 

Take 2, 2;' G C, > 0, < ^ < 1, and define a distribution on the set of all Young 
diagrams by 

M.,.,,,,(A)-(i-oe H^x^e)H'{\ey 

We have used the following notation: t = zz'/6; 

{^)x,o= n {z+{J-l)-{^-m, 

(*j)6A 

where the product is taken over all boxes in a Young diagram A, {i,j) stands for the box 
in ith row and jth column; |A| is the number of boxes in A; 



m,o)= n ((a^-j) + (a;-^)^ + i) 



H\x,e)= n {{x^~J) + {X■-^)e + e), 

where A' denotes the transposed diagram. One can show that J2xeY ~ 
where the sum is over the Y of all Young diagrams. If, for example, z' = z, then all the 
weights are nonnegative, and we obtain a probability distribution on the set of all Young 
diagrams. Mz^z',9,^{X) is called the 2;-measure. Details and explanations of importance of 
^-measures in representation theory can be found in Borodin and Olshanski [6j, Olshanski 

BSl. 



Proposition 4.1. For N = 1,2,... let Y(N) C Y denote the set of diagrams A with 
/(A) < N. Under the bijection between diagrams A G Y(A^) and N-point configurations 
on Z>o defined by 

A < — > XN-i+i = \i-2i + 2N {i = l,...,N) 
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the z-measure with parameters z = 2N, 9 = 2, z' = 2N + /3 — 2 turns into 

Prob {xi, . . . , Xn} = const ■ Y\ p"^^' Y\. ~ ~ ~ ~ + 1); 

i=l l<«<i<A^ 

which is precisely the discrete symplectic ensemble with the Meixner weight in the sense 
of Definition \2.1\ 

Proof. The proof is a straightforward computation based on the apphcation of the exphcit 
formulae for i^(A; 2)H'{X; 2), see the proof of Lemma 3.5 in [6j, and {z)x,e, see Section 1 
in [6]. ' □ 

Proposition 4.2. For N = 1,2, ... let Y'(2N) C Y denote the set of diagrams A with 
/(A') < 2N. Under the bijection between diagrams A G Y'(2A^) and 2N-point configura- 
tions on Z>o defined by 

A ^ X2N-i+i = 2X[-i + 2N {i = l,..., 2N) 

the z-measure with parameters z = —2N, 9 = 2, z' = —2N — (3 turns into 

Prob {a;i, . . . , X2n} = const ■ JJ JJ {xj - Xi), 

where 



i=l l<i<j<2Af 



X. 

and 



2 ■ A ■ . . . ■ Xi, Xi is even, 
1 ■?>■...■ Xi, Xi is odd, 



(xj + /5 — l)(a;j + /5 — 3) . . . (/3 + 1), Xi is even, 
{xi + [3 — l){xi + (3 — ?>) . . . (3 , Xi is odd. 

This is a discrete orthogonal ensemble in the sense of Definition \3. ll 

Proof. The proof is also a straightforward computation based on the formula 

n (2a:-z-2a;.+j) 

1 i<«<j<K^') 



H{X;2)H'{X-2) n!L'?(2A', - ^ + /(A'))! 



□ 



5. The derivation of the correlation kernel for discrete symplectic 

ensembles 

Recall that the Pfaffian of a 2N x 2N antisymmetric matrix A = \\ Aj^ \\'j^=i is defined 



as 

PfA= sgn{a)Ai^i^ . . . A^^^^i^^. 



(T=(ii,...,j2Jv)eS'2JV 
n<«2v:«2JV-l<«2JV 
jl<i3<---<i2iV-l 



One has {Pf Af = det A. 
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Lemma 5.1. Assume that ipi, . . . , ip2N and ipi, . . . , ip2N are arbitrary finitely supported 
functions on Z>o. Set 



L i'2N[ 



and introduce a 2N x 2N antisymmetric matrix A = [Aij].^^^ whose entries, Aij, are 
given by 

We have 



(5.1) 



det [^(xi), V^(xi), . . . ,^(x7v), V^(xAr)] =FiA. 

x={xi<...<xn)C'Z>q 



Proof. This is one of de Bruijn's formulas, see de Bruijn [TTj. 
Lemma 5.2. Let 



iri-Ax) =x' ^ + 



z = l,...,2iV, 



is an arbitrary system of monic polynomials of degrees 0, . . . , 2N — 1. Set 

ipi{x) = 7r,_i(x), iJi{x) = 7ri_i(x + 1). 

Then 

Y[ {Xi - Xjf [{Xi - Xjf - l) = det [^(Xi), V^(Xi), . . . ,^(xAr), V^(xAr)] . 
l<i<j<Af 

Proof. We have 



□ 



l<i<j<N 

l<j<j<Af 

= V{xi + 1, Xi, X2 + 1, X2, . . . , Xat + 1, Xn) 

= {-l)^V{xi, xi + 1, X2, a;2 + 1, . . . , xn, xn + 1' 
= (-l)^(-l)^^^^^ det (7r,_i(a;i), 7r,_i(xi + 1), 
= det [^(xi), V^(xi), . . . ,^(x7v), V^(xAr)] . 



, 7ri_i(xAr),7ri_i(xAr + 1)) 



□ 
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Proof of Proposition 12. 3L 

Applying Lemma 15.11 and Lemma 15.21 we obtain the identity 

= pf [Qij]ij=d , 

where 

(5.2) Qij = w{x) {iTi{x)iTj{x + 1) - iTi{x + l)7rj(a:)) , 

and < i,j < 2N—1. Since the partition function Z^^ of the discrete symplectic ensemble 
under considerations is strictly positive, the determinant of the matrix Q is nonzero, and 
the matrix Q is non-degenerate. The matrix M^^^^ defined by equation fl2.4l) is related to 
Q via Q = diag(||7ri||i/2, . . . , ||7r2;v-i||'/') ■ M ■ diag(||7ri|| ^2, . . . , ||7r2^_i ||i/2). □ 
Proof of Theorem 12.41 

Let ({x) be some finitely supported function defined on Z>o. The following identity holds 

true 

(5.3) 

TV 

{xi<X2<...<xn)CZ>o «=1 ^<i<j<N 

where the matrix elements Aij{() are defined by the formula 

(5.4) Aij{C) = ^ [ni^i{x)nj^i{x + 1) - 7ii^i{x + l)7ij-i{x)]w{x)Cix). 

By the same argument as in Tracy and Widom [?T], §8, we find 



r 2Af-l 2N-1 

E n,ix + l)Q;^\M^'^'i^W^Hy) - E 7r,ix + l)QT^\,iy+l)w'/\x)wy'{y) 

i,j=0 id=0 
27V- 1 2N~1 

E n,{x)Qi/7r,{y)^'/\x)w'/\y) - E Mx)Qij'7,,{y + l)w'/\x)w'/'{y) 

i,j=0 t,j=0 

where the variables x, y take values in Z>o, Qij is defined by equation (15. 2p . and Hi^x) = 
x^ + . . .] z = 0, 1 ... is an arbitrary system of monic polynomials of the discrete variable x. 
It is now straightforward to check that the kernel KN4{x,y) just written above is exactly 
the kernel of the operator K^^, where 

D^S]\i4 —D-^Sn^D^ 
Sna —SnaD- 



K 



Ni 



Finally observe that the matrix Aij{C,) remains unchanged if we replace 7rj_i(a; + 1) by 
7rj_i(x + l)-y^^7rj_i(x), and 7ri_i(x+l) hj -Ki-i^x + l) - ^J^^^^:i-l{x) in equation 
(15. 4p . This results in formula (12.61) . Theorem 12.41 is proved. □ 
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6. The derivation of the correlation kernel for discrete orthogonal 

ensembles 

Lemma 6.1. The probability of a particular configuration xi < . . . < X2n of the discrete 
orthogonal ensemble (see Definition \3. 1\) can be rewritten as 

2N 

(6.1) Prob{Xi, . . . ,X2Ar} = Yl i^i- Xj)'P^H^'h^j)]'i!j=l 

1=1 l<i<j<2N 

where w{x) is defined in terms of the weight function W{x) by formula 1^3. 

Proof. We will compute the Pfaffian in equation fl6.ip . and will show that (16.11) coincides 
with the expression for the probabihty of a particular configuration xi < ... < X2n 
in Definition 13.11 Observe that the semi-infinite matrix e defined by equation (12.31) is 
representable as follows 

e = jrxjr^ 

where 

/(O) 
/(I) 
jr= /(2) 
/(3) 



/(O), /(I), /(2), ... are defined for A; = 0, 1, 2, ... by 
1 w{2)w{A)...w{2k) 



fm 



and 



y^2^w;(l)«;(3)...w;(2A;-l) 



, /(2A;+1) 



w{l)w{3) ...w{2k + l) 



,/w{2k + l) w{2)w{A)...w{2k) 



T 



0-10-10-10-1 

1 

-1 -1 -1 

1 1 
-1 -1 
10 10 10 



That is, T is an antisymmetric matrix whose entries are defined by the relations 
T(2^ + 1, 2j + 1) = T(2z, 2j) = 0, for any z,j > 0, 

T(2i + 1, 2j + 2) = 0, for < i < j, 
T{2i, 2j + 1) = -1, for < i < j. 
Since relation (13. ip implies that 

f "'(l)^'(3)...^W • 11 

W(x)-< «'(2)^(4)...«'(^-l)' ^ISOQQ, 
VV [X) — <. w(2)w(4)...w{x) 

[ w(i)wi3)...w{x-i)^ X IS even, 
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we obtain 
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2N 



]2N 
\i,j=l 



1=1 



J]wi/2(a;,)Pf[e(x„x,)]; 

2N \ / 27V 



'2N 



,j=l 



X{W{xi) Pf 



T(xi, . . . , X2N 



Xl, . . . ,X2Nj 



Xi, . . . , X2Ar) 



The proof is completed by the following Lemma, cf. Borodin and Strahov [9], Section 
3.2. □ 



Proposition 6.2. 



Pf 



Tfxi, . . . , X2N 



Xi, . . . , X2n) 



(—1)^, if Xi — Xi-i is odd, and xi is even, 
0, otherwise. 



Proof. If xi is odd, then the first row of the matrix T(xi, . . . ,X2n 



xi, . . . , X2n) consists 



of only zeros. Thus, if Pf 



T(xi, . . . , X2N 



Xi, . . . , X2n) 



assume that X2i-i and X2i have the same parity. In this case rows 2i — 1 and 2i of 



7^ 0, Xi must be even. Now 



T(xi, . . . , X2N 



Xi, . . . , X2n) are equal to each other. Therefore, if Pf 



Tfxi, . . . , X2N 



Xi, . . . , X2N) 



0, the elements of the set X = (xi, . . . ,X2n) are such that Xi is even, X2 is odd, X3 is 
even, and so on. This proves the condition on the parity for the configurations X = 



(xi, . . . , X2n), for which Pf 



T(xi, . . . , X2N 



Xi, . . . , X2n) 



7^ 0. Moreover, using the defini- 



tion of Pfaffian, it is not hard to conclude that Pf 



T(xi, . . . , X2N 



Xi, . . . , X2n) 



for the configurations for which Pf 



T(xi, . . . , X2N 



Xi, . . . , X2N> 



7^0. 



\N 



□ 



Lemma 6.3. // 0i, . . . , (j)2N o^^e arbitrary finitely supported functions on Z>o, and e(x, y) 
is an antisymmetric function defined on Z>o, then 



(6.2) Yl ^^■)]'j=i = 

a;i<...<a;2jvCZ>o 

Proof. This is another de Bruijn formula, see jllj . 



^ e(x,?/)0j(x)0fc(?/) 



□ 
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Proof of Proposition 13.31 

Consider the partition function Z^i, which is defined by 

2N 

Zm= JJ {xi - Xj)Ff[e{xi,Xj)] 



2N 



xi<...<a;2jv i=l 



l<*<i<2Af 



Write 



2N 



l<^<i<2A^ 



-l)^det [w'/\xk)n,{xk)] 



i=l 



where < j < 2N — 1, 1 < k < 2N, and {ttj^x)} is an arbitrary system of monic 
polynomials. Apply Lemma 16.31 and obtain 

-1 2N-1 



:-l)^det 



^ e(x, y)TTj{x)w^^^{x)TCk{y)w^^'^{y) 

x,y&Z>o 



J k,j=0 



Since Z^i is nonzero, the matrix with the j, k entry 

^ e(x, y)TTj{x)w^/^{x)TTkiy)w^^^iy) 

is non-degenerate. This implies that the matrix M*^^-' whose j, k entry is given by formula 
(13. 2p is non-degenerate as well. □ 

Proof of Theorem 13.41 
Use Lemma 16.11 to rewrite the formula in definition 13.21 as follows 

2N 

Ww'"\^i){^'r^{^^)) n (a^.-a^.) Pf ^otZ=X = V^et (/ + ^Kmi). 

xi<...<X2N «=1 l<«<i<2iV 

The rest of the argument is very similar to the derivation of the correlation kernel in §9 

□ 



of Tracy and Widom [JT], and we omit the details. 

7. The general identities 
Lemma 7.1. We have for < i < 2N — 1 

(7.1) (SNiK^Dipi) (x) = ipiix), {SmKN^^i) (x) = (fi{x 



(7.2) 



S, 



N4 \ni 



\ni 



0. 



where Tij^ denotes the complement ofTii^ inTi. 
Proof, a) For < i < 2N — 1, and x G Z>o we have 



2N-1 



j=0 
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The second sum in the equation above can be rewritten as follows 



2/ez>o 



w{y+l) \j w{y) 

= ^ Pj{y)Pi{y + '^)w{y) - ^ Pj{y)Pi{y - '^)w{y - i) 
= E (pAy)Piiy + ^)-pAy + ^)Pi(y))^(y) 

where {pj} is the family of polynomials ortho normal with respect to the weight w. There- 
fore, 

27V- 1 
j=0 

The action by the operator Sn4 from the left gives 

2Af-l 

(SNiKiyDLfi) (x) = ^ SNiix,y) ^ ipj{y)M, 



2N-1 

(4) 
ji 

>o i=o 



2N-1 2N-1 



E E Vk{x)fj^^ki^<^iiy) E = '^ii^) 

>0 k,l=0 j=0 



b) It is clear from the definition of the operator Sna that (SN^ipj) (x) = for j = 
2N,2N + 1,.... In other words, Sn4 I-^x = 0. 

c) The formulas for Sni can be proved by the same procedure, replacing D everywhere 
by e. □ 

Lemma mi shows that the operators 5'Ar4, Sni, Kn, D and e satisfy the same algebraic 
relations as the corresponding operators in the continuous case considered in Widom |43j . 
Sections 2-4. Therefore Theorem 12.61 and Theorem 13.51 can be established by arguments 
from Widom [13], Sections 2-4. Here we only give a simple "matrix" explanation why 
formulae in Theorem 12.61 and Theorem 13.51 indeed hold true. 

Let us write each operator in the basis yjo, V^i, ■ ■ ■■ Then each operator A which acts 
in Ti has the representation 
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where the upper-left corner corresponds to the subspace TIn spanned by (/jq, • • • , V2N-i- 
In particular, K]<q is representable as 



and we have 



In 




i^)N 0\ _ , ^ ,_^0 (A) 



In what follows we formally manipulate with operators. Equation (17.11) implies 

{Kr,SN,Kr,){Kr,DKN) = Kr,. 
We can rewrite the above identity as follows 

(^7V4)7V \ / {D)n ^ \ = f 

J \ J \ 
Therefore {S]y4)')-i^ = ( ^^g^ g j ? and we obtain 





In order to express {D)q{D)]^^ use the relation 

eD = I. 

This relation can be rewritten as 

(e)o (6)2 ) V P)o {Dh I 

and we obtain {e)N{D)N = In — (e)i(-D)o- Multiplying this equation by {D)q from the 
left we have 

(D)o(e)^(D);v = {D\{In - (e)iP)o) = {In - P)o(e)i)P)o, 

therefore 

P)oP)^^ = (/-p)o(e)i)-^P)o(eV 
We conclude that the following formula holds 

{D)nASN.)n. = ( (/ - {DUe)lY\DUe)^ 
Now we are going to show that 

(7,3) ( _ ) = - f^^^" - '<''OK.yr K.. 
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Indeed, 

DKn - KnDKn = 



I - {DKn - KNDKN)e 





P)o 





and using the formula ^^^^ =^ h^^a equation (17.31) . Thus we 

have shown that 

{D)hASna)h^ = [I- {DKn - KNDKN)e]-' K^. 

Taking into account equation (17. 2p . we obtain the formula in Theorem 12. 6[ The formula 
in Theorem 13.51 can be deduced by similar computations. 

8. Proof of Theorem 12.71 

The idea of the proof is the same as in Widom ^1], section 3. 

Recall that the operator D acts on the elements of T-Cn according to the formula 



N / N / w(x) , , / w(x — 1) , 

Set Ly) = (/ — K]^)D, and agree that the domain of this operator is TIn- Denote by Nd 
the null space of L/), i.e. 

Md = {f\LDf = 0,(f e Hn}- 

We want to find A/"^. 

Let us find the null space J^n of = (/ — K]\f)D. The general element of Ti.^ is of 
the form Lp = pw^^"^, where degp < 2N — 1. Such is an element of A/d if and only if Dip 
is an element of Hn- We have 



V[x + 1 - ^ pix - 1 

w[xj 



It follows that if is an element of Afjj if and only if ^^^^(^^^ pix — 1) is a polynomial of degree 
less or equal to 2N — 1. 

We will denote by ai,a2, ■ ■ - ^ai finite poles of ^^^^^ ■ Let {Ti.N)ai subspace of 

those (f, for which ^l^^^^ pix — 1) = 0(1) in a neighborhood of aj, and let {Hn)oo be the 

subspace of those y?, for which ^l^(^-^^ pix — 1) = 0(x^^~^) as a; — oo. Then 
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and 



= (Ij^^^^i) + ^^^^ 



Let be the finite pole of "^^^^^^^^ of order rij. Observe tliat p{x — 1) "'^^^^^"^-' = 0(1) if and 
only if 

(8.1) p{ai - 1) = 0, p\ai - 1) = 0, ... , p("^-^)(ai - 1) = 0. 

Therefore ip is an element of {7iN)ai if ^"^^ only if condition (18. ip is satisfied. 
Now expand ip in terms of basis elements ipi, . . ., ip2N~i of Hn 

2N-1 
3=0 

Thus 

2N-1 

j=o 

and condition (18. ip implies that (p is an element of {Ti.N)ai ^^"^ only if 

2Ar-l 

(8.2) B,pf\a, -1) = 

j=0 

for < k < rii — 1. Set 

2Ar-l 

Ck= fjpf\ai - 1). 
i=o 

Condition (18.20 implies that G {'HN)ai if and only if is orthogonal to all ^k- Thus, 
^fc's span (7iiv)i- 

Clearly, {TCn)oc is the span of v'a; for < A; < 2A^ — — 1, where rioo is the order of 
'^l^^^^ at oo. Therefore CHn)^ is the span of v^fc for A; > 2N — 71^0 ■ We conclude that the 

I 

dimension of the span of all {Hn)'^. (including {TCn)^) is at most riac + J2^ar Pick an 

i=l 

I 

orthonormal basis ipi, . . . ,ipn of this span. Note that n < n^o + ^ Ua^. 

i=l 

n 

Since -i/^i, . . . , ■i/'n is a basis of A/"^, we have = ^ l^D^i ® ipi- Therefore we obtain 

1=1 



Ld = {I- Kn)D = - KN)Diji ® 



i=l 
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Set ipi = {I — Kjq)D%lJi. Since Li:,tpi 7^ we conclude that iIji G 'Hj^, and we have 

n 

[D, Kn]Kn = ^ ® ipi, i'i e H^, i^i e Hn- 

i=l 

Thus the first part of Theorem 12.71 is proved. 
By the antisymmetry of e we obtain 

n 

1=1 

In order to find D-Hj^Sna we need to compute the inverse of / — [D, Ki^]Ki^e. Using the 
formula 

(/ + ^ a, ® hi)-^ = 1- ^(T-^)i^a, ® bj, 

hj 

where T is the identity matrix plus the matrix of inner products {bi,aj) (assuming T is 
invertible), and ® 6i is a finite rank operator, we obtain the formula for D-h^Sna 
stated in Theorem 12. 7[ □ 



9. Proof of Theorem 13.61 

Set = {I — K]\r)e, and agree that the domain of this operator is T^at. Denote by Af^ 
the null space of L^. Observe that 

AC = {v\v e Hn, et' G T-In}- 

Assume that ip G Af^, and u G AC. Since u G Hn, and Kn is a symmetric operator, 
we have {u,KNeip) = {u,eip). Since u G AC, eu is an element of Hn (see above). But 
Deu = M G H-N, which implies that eu belongs to A/13, see previous Section. Therefore, 
{u,eip) = —{eu,ip) = 0, and 

[u, K^etp) = 0, for all G Af^ and u G A^. 

We conclude that Kj^eKjsi takes into A^-^. 

Let us show that A^^ has the same dimension as M^. The operators Kj^eKjsi and 
KpfDK^ are invertible because we have already shown that their matrices Mj)^ and M^'^ 
in the basis ipo, . . . , ip2N-i are invertible, see Propositions 12.31 and 13. 3[ We just saw that 
K]yeK]\f takes Af^ to A^^, and very similar arguments show that Kp^DK]^ takes A^"^ to 
M^. Therefore dimA/"^ = dimA^^. 

Now take ipi, . . . jipn found in the previous section and orthonormalize the functions 
K^eipj. As a result obtain an orthonormal basis in A^^. Denote this basis r/i, . . . , 

This gives 

n 

L,= {I- KN)e = ^(/ - KN)er]i ® 77^. 
1=1 
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Set Tji = {I — Ki\i)erii. By the same argument as in the previous section we find (assuming 
invertibihty of U) that 

n 

e-HMSm = Kn- ^{U^^)ij{fli) ® {KNDrij) , 

where the matrix U is defined by its matrix elements by 

Uij = 5ij + {Drji^fjj). 

This immediately gives the formula in the statement of Theorem I3.6[ □ 

10. Proofs of Corollary 12.81 and Corollary 13.71 

Proposition 10.1. // the commutation relation between the operators D and Kjs[ takes 
the form as in the statement of Corollary \2.8[ with ipi G 71^, and ip2 £ '^nj then eipi is 
an element ofHN- 

Proof. The explicit form of [D,Kn] in the statement of Corollary 12.81 implies 
(10.1) [e, Kn] = Keipi ® eip2 + eip2 ® eipi). 

(We have used the fact that eD = 1, and the antisymmetry of e). Acting by [e, i^iv] on 
ipi we find 

This gives 

KNeipi = (1 + X{ip2,eipi))etpi. 
Therefore, either eipi is in Hn, or K^eipi = 0. To rule out the second possibility observe 
that [D, Kn]Kn = A?/'2®V'i) it follows from the expression for [D, Kn] in the hypothesis 
of Corollary 12.81 Here 1^2 G Hj^, and ipi G Hn, and Theorem 12.71 implies that ipi is an 
element of Afj^. It was shown in the proof of Theorem 13.61 that K^^ bijectively maps Afj^ 
into M^. Therefore, K^eipi must be an element of Af^, so K^eipi cannot be zero. □ 

Proof of Corollary 12.81 If the conditions of Corollary 12.81 are satisfied we have 

[D, KN\KNe = -Xip2 ® etpi. 

Then Theorem 12.61 implies 

D-HmSna = {I-Hn+dhn + ^'^2 ® eiJiy^K^. 
Now we use the formula 

{I + a^by^ = I - T'^a®h, 

where a = \ip2, b = eipi, and T = 1 + X{ip2, eV'i)- By Proposition 110. II eih-i E Hn, and we 
conclude that T = 1. Therefore 

or 

Sna = ^Kn — Xeil)2 ® eipi. 
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□ 

Proof of Corollary 13.71 

By Theorem 13.61 we need to compute inverse of 

I-[e,KM]K^D. 

From equation (110. 11) we find 

[e, Kn]KnD = -A(e^i DKNetp2 + e^2 ® DKNetpi). 

Now we compute DK^eipi as follows 

DK^e^Pi = [D, KN]eiJi + K^Dei^i = [D, fC^JeV'i + 

= A(V^i (S) ip2 + ip2 (S) ^i)eipi + 1^1 = i>i. 

(We have used {ipi^eipi) = 0, and the fact that by Proposition 110.11 {ip2,eipi) = 0). In a 
similar way we compute DKNeip2 

DKr,eij2 = [D, iT^je^'s + K^Deij2 = [D, Kr,]eij2 

Therefore we have obtained 

/ - [e, Kn]KnD = I + \etp2 ® ^i- 

Using the formula 

{I + a0b)-^ = I -T~^a0b 

where now a = \e1p2, b = ipi, 

T=l + {b,a) = l + {iJi, \eiJ2) = 1 - A(e7Ai, ^2) = 1, 

we obtain 

(J - [e, Kn]KnD)-^ = I- \ei^2 ® V-i- 

This gives 

^Hn^ni = Kn - Ae^2 ® V^i, 
and we finally arrive to the formula 

□ 

11. Discrete Riemann-Hilbert Problems (DRHP) and difference 
equations for orthogonal polynomials. 

In this section we recall a few claims from Borodin and Boyarchenko [5]. Let X be a 
discrete locally finite subset of C Let w : X — >■ C be a function. Assume that all moments 
of w are finite. Denote by the space of polynomials in the complex variable (, and 
introduce the inner product in C[(] 

(/(C),^7(C)L:=E/(^)^(^M^)- 
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If the restriction of (.,.)«) to the space CfC]-'^ of the polynomials of degree at most d is 
non- degenerate for all ci > 0, then there exists a unique collection of monic orthogonal 
polynomials {P„(C)}^o associated to w such that if m 7^ n, then {Pn{C)i Pm{C))w = 0, 
and {PniQi Pn{C})w 7^ for all n. If this condition holds, we say that the weight function 
w is non-degenerate. Note that if X is finite, and consists of + 1 points, the inner 
product (., is necessarily degenerate on C[C]-'^ for all d > N. In this case we require 
that (., .)w be non-degenerate on C[C]-'^ for < d < N, and we are only interested in a 
collection of orthogonal poynomials of degrees up to N. 
As in Refs. [SI IH [5] , we say that an analytic function 

m : C\X ^ Mat{2, C) 

solves the DRHP (X, w) if m has simple poles at the points of X and its residues at these 
points are given by the jump (or residue) condition 

Res m{() = hm {m{()zu{x)) , x G X, 

where 

, . f w{x) \ 

Theorem 11.1. Let {Pn{C)}n=o collection of monic orthogonal polynomials associ- 

ated to w, where N = card{X) — 1 G Z>oU {00}. For any k = 1,2, . . . , N the DRHP(X, w) 
has a unique solution mx{C) satisfying the asymptotic condition 

where I is the identity matrix. If we write 

^^^^^ - ^ m|(C) mf (0 ) ' 

thenm'i{C) = PkiC), m|(C) = (P.-i, P.-i).^Pfc-i(C), m|(C) = E andmf (() = 

xex 

-rfc-ij^ 2^ . 

x€X 

Proof. See Borodin and Boyarchenko [5], Section 2.3. □ 

Let X be a finite or a locally finite subset of M, card(X) = + 1, where A^ G Z>o U 
{00}. Assume that X is parameterized as X = {'n'x}x=o^ ^^"^ that there exists an affine 
transformation cr : C — C such that cttix+i = tTx for all < x < A^. Denote the derivative 
of a (which is constant) by rj, so that 

a(Ci) - a(C2) = viCi - C2) for all Ci, C2 G C. 
Let w be a strictly positive real valued function defined on X. Then w is non-degenerate. 
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Theorem 11.2. Assume that there exist entire functions di{C) and d2{C) such that 

«;(7r^_i) di{7T,^) 

di{iTo) = 0, 
d2{(J~^H]y) = 0, if N is finite. 

a) We have 

wherec= {Pk-i, Pk-i)7j , and M^^ {() , M^"^ {() , M"^^ {() , M"^^ {() are entire functions. More- 
over, the matrix M{() whose entries are M^^{(), M^^{(), M'^^{(), M^^(C) is given by 

where 

■ rfi(C) 
d^iO 



and mx{C) is the solution of the DRHP(X,w) with the asymptotic condition 

b) If it is known that di{() and c?2(C) ore polynomials of degree at most n in (, and 

di{C) = -^iC" + {lower terms), 

d2{() = A2C" + {lower terms), 
then 2x2 matrix M{() is a polynomial of degree at most n in (, and 

M'HC) M^'{0\_(v'Xi + (lower terms) 

M^\C) M^\C))~\ ^-k^^ jQ +{^ower terms). 

Proof. For the proof of these statements see Borodin and Boyarchenko [5], Sections 3.1- 
3.3. □ 

12. Commutation relations 

In this section set X = Z>o, and crC = C ~ 1 fo^' ^'^Y C ^ Assume that w{x) is 
a strictly positive function on Z>o, which satisfies the condition in Theorem 111.21 i.e. 
the ratio of w(x — 1) and w{x) equals the ratio of entire functions di{x) and ^2(2;), and 
cii(O) = 0. 



DISCRETE ENSEMBLES 27 

Proposition 12.1. We have for n = 1,2, . . . 

[D_LPn) [X) = Lpn{x) H -75 -j:r(^n~l{x) , 

d2{x){Pn,Pn)'J\Pn-l,Pn-l)'J' 

[D_^n-l) {X) = — ipn{x) + —-—fn-l{x). 

d2{X) d2{X) 

Proof. These equations are equivalent to the difference equations for the monic orthogonal 
polynomials in Theorem II 1.21 a). □ 

Proposition 12.2. Introduce the 2x2 matrices D±{x) by the formula 



D±if2N^l{x) J \ Df{x) Df{x) J \ Lp2N~l{x) 

The matrix elements of D^{x) are related with the matrix elements of D_{x + 1) as 
:i2.2) 



D\_\x) = Di\x + 1), D'^{x) = -Dt'{x + 1) 



^21 (3,) ^ -Dl\x + 1), Df{x) = Dl\x + 1) 
Proof. The definition of D± implies that 

(12.3) D+{x)D4x + l] - "^^^^ 



W[X + 1) 



Observe that detM(x) = di{x)-d2{x) (see [5], §3). This fact together with the assumption 
that the ratio of w{x — 1) and w{x) equals the ratio of di{x) and d2{x) imply that the 
determinant of the matrix D^{x + 1) equals which (together with equation (112.31) ) 

leads to the relation between the matrix elements in the statement of the Proposition. □ 

Proposition 12.3. We have 
[D,K^]{x,y) 



Dl^{x+l)-Dl^iy) D22(x+1)-D?2(y) 
x+l—y x+l—y 



a2N{V2N{x),ip2N~l{x)) I gnffl)^D"fa) DlH^+^l)-DiHy) I ( y?t^^|^y) 

\ x+l-y x+l-y / ^ 

/ Dl\x)-Dl\y+1) Dl\x)-D]}(y+l) \ / ^ ^y^^ 

+ a2N{V2N{x),^2N-l{x)) \ Di2(4:|rjfa+i) DlHx)~-D^} {y+1) 1 ( ^^^_^{y) 

\ x—y—1 x—y—1 / 

where a2N is the coefficient in the Christoffel-Darboux formula for the kernel K{x,y), 

K [x, y) = a2N • 

x-y 
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Proof. A direct computation of the commutation relation between the operators D and 
K gives 

[D,K\{x,y) = a2N 

I X + 1 — y 

{D_ip2N-l) {y)^P2N{x) - {D_(P2n) {y)'f2N-l{x) 





X + 1 


-y 


{D_LP2n) 


{x)^2N-i{y) - 


{D_ip2N-l) {x)(f2N{y) 




X — 1 


- y 


{D+(P2N- 


i) {y)(p2N{x) - 


{D+if2N) {y)ip2N-i{,x) 




X — 1 


-y 



+ 

Express D±(p2N, D±f2N-i as hnear combinations of (f2N and ip2N-i whose coefficients are 
entries of matrices -D+(x), In the resuhing formula replace the matrix elements 

D+(a;), D^^{x), Dl\x), Df{x) by the matrix elements Dl\x + 1), D}_^{x + 1), D^^x + l), 
D'^{x + 1) in accordance with formula fll2.2p . This gives the formula in the statement of 
the Proposition. □ 

Proposition 12.4. Assume that the functions di, d2 are polynomials of degree at most 
m. Furthermore, assume that d2 has zeros at points ai, . . . ,ai of degrees ni, . . . ,ni corre- 

spondmgly. ihen each of the expressions x+i-y ' x+i-y ' x+i-y ' 

D]_^{x+l)-D^J{y) ■ n ■, r ■ f+U f 

^j^^_y IS a finite sum of expressions of the form 

^m—l—n\ — ...—ni I rii „ \ /m— 1— ni — ri; I Ui „ 



A:=0 i=l k , 

and each of the expressions 



„ _ „ ... — n 

k=Q i=l fc,=l ^ ' / V k=Q i 



^ (y — a,)* 
1 ki=\ ^ 

D2l(x)-Dil{j/+1) D^}(x)-D^}{y^\) Dl^{x)-Dl^(y+1) Dl^(x)-Dl^{y+1) 



x—y—1 ' x—y—1 ' x—y—1 ' x~y—l 

is a finite sum of expressions of the form 

^m— 1— ni — n; I Ui ^ \ /m— 1— ni — n; I rii ^ 

E ■^»-' + EE7r^ E ^'■/+EEfeTr^ 

fc=0 i=l k,=l ^ ^' J \ k=0 i=l ki=l 

where Ak, Ck, Bk^ and Dk^, Ak, Ck, Bk^ and Dk^ are some constant coefficients. 

Proof. If di{x) and ^2(2;) are polynomials of degree at most m, then M^^(x), M^^(x), 
M^^(a;), and M^^(a;) are polynomials of degree at most m, see Theorem \11.2\ b). By 
Proposition I12.ll and by the very definition of the matrix D_{x) each D^}{x), D'^{x), 
D}^{x), D'^{x) is a polynomial of degree at most m divided by d2{x). Set 

d2{x) = const -{x — Oi)"^ . . . (x — ct;)"', 

where 

ni + n2 + ■ ■ ■ + ni < m. 
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Denote D}}{x)d2{x) = AmX"^ + . . . + Aq. Then we have 

Di\x + 1) - Di\y) ^ Amix + l)"' + ... + Ai{x + l) + Ao 
x + l-y ~ (x + 1 - (x + 1 - ai)"i . . . (a; + 1 - a;)"' 

ArnV"" + ... + Aiy + Ao 



Since 



ix + l-y)iy-a,)-^...iy-ai)^i' 

{x + 1)^ - y^ 



X + 1 — y 

^ a+b<k~l 

we arrive at the result. The same considerations are apphcable to all expressions in the 
statement of the Proposition. □ 

13. Difference equations for the orthonormal functions associated to 

THE MEIXNER and TO THE ChARLIER WEIGHTS 

Proposition 13.1. Set 

where WMeixner{x) is the Meixner weight defined by equation ^2. 9\) . and {Pj{x)} is the 
family of the monic Meixner polynomials. The functions {^ni.x)}'^^^ satisfy the following 
system of difference equations 



^JCX[X + /? - \)^n[x - 1) = (X - n)ifn{x) + Cuin + /? - \)^n-\[x), 



y/cx{x + (3 - l)(fn-i{x - 1) = c{x + {u + (3 - l))(fn~i{x) - cn{n + (3 - l)(fn{x). 

Proof. Let m{() denote the solution of DRB.P{Z>q,w Meixner) satisfying the asymptotic 
condition 



see Section [TTJ The matrix m{() has a full asymptotic expansion in ^ as C — * oo, and we 
can write 

(13.1) ^(0(V ])c' + o{C'). 

By Theorem 111.21 there exists an entire matrix- valued function M{() such that 

(13.2) M{0=m{C-l)D{Om-\0, 
where the matrix D[() has the following form 



C 

cC + ci(3-l) 
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We can rewrite equation (113.21) as 



M(C) 



m(C - 1) 



(c-1)- 








(C-1) 



C 

cC + c(/3-l) 



(1 

The last multiplier in the expression above has the form 



m(C) 



C^" 
C" 



-1 



+ 0{C 



C-n 




0(C 



1 - aC^ 

the multiplication of the two matrices in the middle gives 

cC + cn + c{(3 - I) ) 

and the result of the multiplication of the first two matrices in the expression for M{Q is 
as it is evident from (113. ip . The straightforward calculation gives 

C-^ (c-i)r VorrM 

[c-l)7 cC + c{n + p-l) J^^^"^ >■ 

Since M((^) is entire the last term 0(C~^) is identically zero by Liouville's theorem. By 
Theorem 1 1 1 . 2 1 the monic Meixner polynomials satisfy the following system of the difference 
equations 

CP„(C - 1) = (C - ^)Pn(C) + c„-i(c - l)rP„_i(C) 

C„-lCPn-l(C - 1) = -(C - l)7Pn(C) + C„-l (cC + c(n + /5 - 1)) P„_i(C) 



M(C) 



(13.3) 



where c„ = (P„,P„)^ . Moreover, we also have the condition detM(^) = detD(C) which 
follows from the fact that the determinant of the solution of the Discrete Riemann-Hilbert 
problem identically equals 1, and from the relation between M(C) and D{(), equation 
(113.21) . This condition implies a relation between parameters, 

(13.4) (c-l)V7 = cn(n + /?-l). 

Replacing n — 1 by n in the second equation of system (113. 3p . and inserting the result into 
the first equation we obtain the recurrence relation for the monic Meixner polynomials. 
Since the recurrence equation for the Meixner polynomials is known (see, for example, 
equation (1.9.3) in Ref. ^25j) this (together with formula (113.40 ) determines coefficients 7 
and r in system (113. 3p . Finally, rewriting (113.30 in terms of functions (pn{x) and ipn-i{x) 
we obtain the difference equations in the statement of the Proposition. □ 
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Proposition 13.2. Let {fn{x)}^^Q be the family of the orthonormal functions associated 
with the Charlier weight defined by equation \2.1(A The functions {v^n(a^)}^o satisfy the 
following system of difference equations 

yJaxLpJyx - 1) = (x - n)ipn{x) + ^-A^) , 
-JaxLpn-\{^ - 1) = -Van^Pnix) + a(fn-i{x). 

Proof. The system of difference equations for the orthonormal functions {v9„(x)}^o asso- 
ciated to the Charlier weight follows from the corresponding system of difference equations 
for the case of the Meixner weight. Indeed, we have the following limiting relation 

(13.5) lim M„ (x; P, = C„(x; a) 

between the nth Charlier polynomial Cn(x; a) and the nth Meixner polynomial Mn{x\ (3, c), 
see, for example, Ref. |23j. □ 



14. The Meixner and Charlier symplectic and orthogonal ensembles 



By Propositions 112.31 112.41 the kernel of the operator [D, Kn] is expressible in terms of 
the following functions: 



(14.1) 



x^'^2N-i{x), x^Lp2N'-, 0<A;<m-l-^ni, 



i=l 



and, for each zero of ^2(0;), the functions 

(14.2) {x + l-ai)-^'^2N-i{x), {x + l-ai)~^^^2N{x)] 1 < ki < rii, 

(14.3) (x - ai)~''*(p2N-i{x), {x - ai)~^V2Af(a;); I < h < rii. 
Proposition 14.1. In the case of the Meixner weight we have 

[D,KM\[x,y) = ^— ^ {C2N{x),C2N-l{x),r]2N{x),r]2N-l{x)) 

( ^2Af(27V+/3-l) -2N^c \ / C2N{y) \ 

y/2N{2N+l3-l) 





(14.4) 







V2N{y) 

\ ri2N-i{y) I 



\ -2N^ y/2N{2N+^-l) 

where we have introduced Cj{x) = and rij{x) = ■ 

Proof. We use Proposition 112.31 which determines the commutation relation between the 
operators D and Kj^i in terms of the matrix elements of the matrix D_{x) defined by 
equation (112. ip . This matrix can be written explicitly in the case of the Meixner weight 
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since we have obtained in Proposition 113.11 the system of difference equations for the 
orthonormal functions '^2N-i{x) and '^2n{.x). Namely, 



D_{x) 



y/2Ni2N+l3-l) 

v^(x+/3-l) 
y/2N{2N+f3-l) x+{2N+t3-l) 
v^{a;+/3-l) x+f3-l 



X-2N 
c(x+l3-l) 



Inserting the matrix elements of the matrix D_{x) into the formula in the statement of 
Proposition 112.31 and taking into account that the coefficient a2N corresponding to the 

case of the Meixner weight equals ^) obtain the formula for [D, K^]. □ 

We want to construct linear combinations ipi,4'2,4'3, i^i of (2N, (2N-i,V2N, V2N-1 in such 
a way that ipi and ^'2 will be lying in and '?/'3, ip^ will be lying in Hj^. 

Proposition 14.2. Set 

(14.5) Vi = V2dVC2N - V/3 + 27V-1C27V-1, 



(14.6) V2 = V2NFi-2N + 1, /3 - 1; /?; c)rj2N - ^cip + 2N - l)F{-2N, p - 1; (3; c)mN^u 



(14.7) 

V'a = (/3+27V)F(/3, /5+2iV-l; [3+2N; c)C2N -^/2cN{p + 2N - 1)F{(3, p+2N; (3+2N+1; c)C; 



(14., 



ipi = ^JP + 2N - lr]2N - V2^mN-i- 



Then ipi,ip2,ip'i,'>Pi are linear independent, ipi^ip2 are elements ofTiN, and ipz^ipi are 
elements of H^. 

Proof, a) Suppose that the linear combination 

,<f2N{x) , (p2N-l{x) 



+ B- 



X + P X + P 

is an element of 7^ at. Here A, B are some coefficients. If ipi is lying in ?-^jv, then 
(14.9) Ap2n{-P) + Bp2N-i{-P) = 

Therefore, we need to check that the equality 

P2iv-i(-/?) 



;i4.io) 



P2n{-P) 



2cN 
P + 2N-1 



is satisfied in order to conclude that ipi defined by (114. 5p is lying in Hn- To check (114.101) 
note that 

Mn{x;p,c) 



Pn[X) 



\Mn{.;P,c) 
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where Mn{x; P, c) is the Meixner polynomiaL We have Mn{x; /?, c) = F(—n, —x; (3] — ) 
see, for example, Refs. [23l ESj for the basic properties of the Meixner polynomials. The 
norm, ||Mn,(.; /3, c)||, is equal to 



Taking into account the formula 

F{-n,p-p--z) = {l + z) 



n 



we see that relation fll4.10p indeed holds. 

b) In the same way we show that '?/'2 defined by equation (114.61) is an element of I-Ln- 

c) Now we need to prove that the functions ipsyip^ defined by equations (114.71) . (114.81) 
are elements of Ti.^ . To show this we note that the space Hat in the case of the Meixner 
weight can be understood as that spanned by the functions 

V^o; ipi] {X + P){X + P - l)x''^yWMeixner{x), < k < 2N - 3. 

All linear combinations of ^^f^^^, ^|^^) '^l'^^^}^^ x+p'-i certainly orthogonal to 

{x + (3){x + (3 — l)x^^y w Meixner (x) , < < 2N — 3. Therefore, we can take as ips^ipA 
linear combinations of the form 

(p2N-lix) . (p2N^ix) ^2Nix) 

^ 7, r ^ TT + ^ 7, 7 + ^ 7, 7) 

X + f3 X + (3 x + (3-l x + (3-l 

where the coefficients C, E and F are subjected to the conditions 
(14.11) 

C -^2^^) + D (^„, "ElM) + E ^^^) + F •^-(^) 



x + P J V x + pj \ x + P-lJ \ x + p- 

(14.12) 

C ^211^] + D ^] + E ^^!^) + F 1 = 0. 



x + P J \ x + P J \ x + P — lJ \ x + P 
In particular, we can take 

i^3 = C^ + D^, ij, = E +F ^'"^ 



x + P x + P x + p-l x + p-l 

provided that the following conditions on the coefficients C,D,E and F are satisfied 
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(It is not hard to check that for n > 1 
(14.14) 



X + 13 



- ^0, -n 

C \ X + p 



C + /3 



x + p 



Taking into account these relations we see that if conditions (114. 13p are satisfied, equations 
ffim]) . (fm2D hold as well). 

Thus it remains to check that for ri > 1 the following relations are valid: 



;i4.i5) 



) _ v^m(/? + n- l)F{p, P + n]P + n + l]c) 



y>,i-i 



{p + n)F{p, P + n - I- P + n]c) 



;i4.16) 



Let us compute the scalar product in the left-hand side of equation (114.151) . We have 

+00 



x + p 



.^_^^pJ{P)n^^M^{x;P,c) iPU^ 



n\ x + P 



X\ 



It is convenient to exploit the discrete Rodrigues formula for the Meixner polynomials 
(see, for example, Ismail [23], Section 6.1): 



Mn{x;p,c) 



{P + n),c^ 



xl 



where 



This gives 



(V/) (x) = /(x)-/(x-l). 



:i4.17) y.0 



'1-c 



,/3 



X + P^ 

where we have used the formula 



n\ ^ — ' X + P 



{x-k)\ 



(vv) (x) = g(^)(-i)V(^-fc). 



Two sums in the righthand side of equation (114. 17p can be rewritten further as 



+ 00 



1 iP + n),c- 



^ k I ^ X + p + k x\ 
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Taking into account the formula 

n 

E(-i 



fc=0 



k I n 
k 



X + k x{x + 1) . . . (x + n) 



we sum up over k and obtain 



' X + (3 



'l-c] 



E 

x=0 



{x + (3){x + (3 + l)...{x + (3 + n) x\ 



:i - of 

/3 + n 



c-^\I^W^FiP,P + n;P + n + l;c). 



T{(3 + n) 



This formula implies that relation (114. ISp indeed holds. Relation (114.161) can be checked 
in the same way. 

□ 

Proposition 14.3. In the case of the Meixner weight the commutation relation between 
D and Kn can he expressed as 



:i4.i8) 



^2N{2N + ,,^01 



^1 

^"4 



where ipi E Hn and ip^ G Hj^ are defined by equations jlj-^ , jl-j-^ correspondingly. 
Proof. Equations (I14.5I) -( |14.8I) can be rewritten as 

/ mil "^12 \ / C2N \ 



;i4.19) 



where 







/ 1 


^2 





















m2i 

77133 "^34 

y m43 m44 J 



C2N-1 
r]2N 

\ V2N-1 ) 



mu = v^ciV, mi2 = -^/p + 2N -1; 
m2i = {p + 2N)F{P, P + 2N -l■p + 2N■c)■ 



m33 



77122 = -^2cN{l3 + 2N - 1)F{P, p + 2N;p + 2N + l;cy, 
^F{~2N + 1,(3-1- (3- c), m^A = c{(3 + 2N - l)F{-2N, /3 - 1; /?; c); 



m43 = y^^+^iV^, m44 = -\/2cN. 
Introduce the following notation 

Ml 



B 



( mn 


mi2 \ 


, M2 = 


' rn33 




V ^21 


m22 J 




\ m43 


TUu J 



^2N{2N + (3-1) 
c- 1 



hi bi2 

^21 ^22 
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2N-1+I3 



622 = ^/2N{2N + P-l), 



where bu = ^/2N{2N + (3 - 1), bu = -2N^, 621 = - 
and 

/ 1 \ 
10 
10 

\ 1 y 

Using equation (114. 19p we express (2N, C2N-1, V2N and rj2N-i in terms of 4'i,4'2, "^3 and ip4, 
and rewrite formula (114.41) as follows 



Q 



Now we compute the matrix {M^)~^BM2^ explicitly. We have 
^1 _ ^2N{2N + (3-l) f 

77122 -^21 



Q 



V^4 / 



(c-l)AiA, 



-mi2 mil 



bu bi2 
^21 ^22 



77144 "^^34 
-m43 77133 



where Ai = det Mi, A2 = det M2. We observe that 



6ll7?l44 


- &12"^43 


= 0, 


&22"^43 


- &2l"^44 


= 0, 


6ll7?li2 


- b2imu 


= 0, 


622^11 


- bi2mi2 


= 0. 



Taking the relations above into account we find 



-1 _ ^/2N{2N + (3 -1) f -(77122&11 - m2ib2i)mu + (m22&i2 - "^2i&: 



(M/ Y'BMo , , . , 

^ ' (c-l)AiA2 

Now we will show that 



-(77I22&II - m2l62l)"^34 + ("^22&12 " "^21&22)"^33 

Indeed, the straightforward algebra gives 

- ("^22&ll - "^21&2l)"^34 + (77722&12 " "i21&22)"^33 



•22)ms3 





A1A2 



= ^/P + 2N-1 [{p + 2N)F{P, P + 2N-l;P + 2N;c)- 2cNF{(3, {3 + 2N- {3 + 2N + I- c)] 
X [{13 + 2N - l)F(-2iV, /5 - 1; /3; c) - 2NF{-2N + 1, (3 - 1; (3; c)] 
From the other hand, 



Ai = ^/p + 2N-l [{p + 2N)F{P, P + 2N-l;P + 2N;c)- 2cNF{P, P + 2N; P + 2N + I; c)] 
A2 = [iP + 2N- 1)F{-2N, p-l-p-c)- 2NF{-2N + 1, /3 - 1; /3; c)] . 
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Thus 

. _ ^2N{2N + /3-l) ( 1 



Formula (114.181) follows after simple algebra. □ 

Proof of Theorem 12. 9|, a) and Theorem 13.81 a). 

Formulae in the statements of the Theorem l2.9[ a) and Theorem l3.8l a) follow immediately 
from Corollary 12.81 Corollary 13. 7[ and Proposition 114.31 □ 

Proposition 14.4. In the case of the Charlier weight defined by equation 1^2. 1C\) we have 



(14.20) [D, Kn] = Y {V2N~1 ® V2N + <^2N (S) ip2N-l) ■ 

Proof. According to Proposition (112. 3p the commutation relation between the operators 
D and Kn is determined by matrix elements of 2 x 2 matrix Z}_ defined by equation 
(112.11) . This matrix can be found explicitly from the system of difference equations for 
the orthonormal functions {<^n{x)}'^=Q associated with the Charlier weight obtained in 
Proposition 113.21 The result is 



Djx) 



X-2N V2N 



1 



Inserting the matrix elements of the matrix D-[x) into the formula in the statement of 
Proposition (112. 3p . and taking into account that the coefficient in this formula equals 
— v2iVa we obtain the desired result. □ 

Proof of Theorem [Ml b) and Theorem [SH], b). 

Since (p2N-i ^ 'Hn, and (p2N £ 'Hj^ we can directly apply Corollary 12.81 and Corollary 
1X71 □ 

15. A LIMITING RELATION BETWEEN MEIXNER SYMPLECTIC AND ORTHOGONAL 
ENSEMBLES, AND THE ChARLIER SYMPLECTIC AND ORTHOGONAL ENSEMBLES 



Theorem 15.1. As l3 ^ oo and c = the correlation kernels for the Meixner symplec- 

tic and orthogonal ensembles with weight -^^c^ (given in Theorem \2.^A a), Theorem \3.8[. 
a) respectively) turn into the correlation kernels for the Charlier symplectic and orthogonal 
ensembles with weight ^ (given in in Theorem \2.9[ b), Theorem \3.S\ b) respectively). 

Proof. In order to prove Theorem 115.11 we apply the formula (113. 5p . If yj^^s*^"^'' denotes 
the nth orthonormal function associated with the Meixner polynomials, and if 
denotes the nth orthonormal function associated with the Charlier polynomials formula 
(113. 5p implies 

hm ™(a;; (3, -^) = ^^'^-'-''(x; a). 
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Therefore if c = -g^, and /5 — > cxd, 

Now the second statement of the Theorem can be checked immediately, taking into ac- 
count that the kernels of the operators e and D in the case of the Meixner ensemble 
become indistinguishable from the kernels of the corresponding operators in the case of 
the Charlier ensemble, as /3 — cx) and c = □ 

16. a limiting relation between the gorrelation fungtions of the 
Meixner and the Laguerre symplectic ensembles 

Fix the measure aLaguerre on R>q defined by 

Consider the set Conf7v(IR>o) consisting of A^-point configurations X = (xi, . . . ,xn)- On 
this set we define a probability measure P^^^ as follows 

Pi^lidX) = const^4 \V{X)\'at^^^UdX), 

where af^g^^^^^{dX) = U^^^OiLaguerreidxi), V{X) = Ui<i<j<Ni^i - ^j) is the Vander- 
monde determinant, and const 7V4 is the normalization constant. 

If we interpret the points xi,X2, . . . ,xn of the random point configuration X as the 
eigenvalues of a iV x quaternion real matrix then the measure Pj^l determines the 
symplectic Laguerre ensemble of Random Matrix Theory, see, for example, Mehta |26j . 
Forrester [T7] . 

Let {Li"'*}^Q be the family of the Laguerre polynomials defined by the orthogonality 
relation 

oo 

J x-e-^Lt\x)L^:\x)dx = + ; + 



and set 



(16.1) ^^:\x) = J 4")(x)xte-t. 

y 1 (a + n + 1) 

For the symplectic Laguerre ensemble the correlation kernel is the kernel of the operator 
K^^l, which is of the form (see, for example, Widom [53]) 
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Here the operator T>S^^\ has the kernel 



(16.3) = + V2^(2^ + ") (V2iVT^C^^^(x) - y^dli(^)) 

and the kernels of S'^^, S^^lv, and VS^^lv can be obtained by action of T> and In 
the formulae written above the function K^^^^x, y) (which is the correlation kernel for the 
unitary Laguerre ensemble of Random Matrix Theory) is given by the formula 

X y 

T) is the operator of the differentiation, £ is the operator with the kernel S{x,y) = 

|sgn(x — y), and Ck"\^) ~ ^~^Vk'\^)- Note that our notation is slightly different from 
that of Widom [43j. 

It is clear from the definition of the symplectic Laguerre ensemble, and from the defi- 
nition of the symplectic Meixner ensemble that there is a limiting relation between cor- 
relation functions of these ensembles. Namely, denote by p{xi, . . . ,Xm) the correlation 
function of the Meixner symplectic ensemble with weight WMeixner{x) given by formula 
12.91 and by pm\xi, . . . , Xm) the correlation function of the Laguerre symplectic ensemble 
defined by weight x"e~^, x > 0. Set P = 1 + a. Then 

(16.4) \im j^^^^Pn. (yT^, . . . , ^) = P^^ (xi, ...,Xm). 

Our aim here is to check equation (116.41) on the level of the correlation kernels. 
Theorem 16.1. Take f3 = a + 1. Then for any strictly positive integers x,y we have 

lim -^DSN^i—^,-^) =VS^^l{x,y), 

c^l- 1 — C V1~C 1 — C/ 



(y^-yt^) =^I'S!?4'(!/.^) = ^S£"4'I'(^.!/). 



40 ALEXEI BORODIN AND EUGENE STRAHOV 

Sketch of the proof. The Meixner polynomials 

M„(x; (3, c) = F{-n, -x; /?; 1 - c'^) 

are related with the Laguerre polynomials L^n\x) via the formula 

hm M„ f-^; a + 1, c) = 
In addition, noting that 



as c ^ 1, 



r(i + ^) r(/3)r(i + ^) Vi-^y r(/3) 

it is not hard to check that 

(16.5) iini^^<^„(^^;a + l,c) = ^^:\x), 

where is defined by formula fll6.ip . Now observe that J il)^^\x)dx = 0, i.e. 



y +00 y +00 



y y 



Indeed, for a > i^'i/'i is an element of iJ^N spanned by v^q , V^i , • ■ • , V^2Af-i)' 

+00 

and we conclude that 8il)^i\'S) = for a > 0. But £il)'f\Qi) = — | / (x)(ia;. To check 



'^^ (a) 

that / ipi {x)dx = holds true for —1 < a < we can use the analytic continuation 



with respect to the parameter a. Using these formulae we arrive to the following 
Lemma 16.2. For any two positive integers x, y 

}}^_Y^Kn ^) = K^j^kx,y), 



c- 
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Using the limiting relations in the Lemma just stated above, and the expression for 
DSn4 in Theorem 12.91 a) we obtain 

lim T^DSN^i—^,-^) = VS^^l{x,y). 

c^l- 1 — C \1 — C 1 — C J 

This is the first limiting relation between the kernels of the Meixner symplectic ensemble 
and the Laguerre symplectic ensemble stated in the Theorem. Other limiting relations 
can be obtained by action of operators V-t and e, and by application of Lemma 116. 2[ □ 



17. Correlation functions for the Meixner orthogonal ensemble and the 

PARITY RESPECTING CORRELATIONS FOR THE LAGUERRE ORTHOGONAL 

ENSEMBLE 

Consider the Laguerre orthogonal ensemble. This ensemble can be defined by the 
probability density function 



2N 



const ■ Y\ 



, 2 



i=l l<j<k<2N 



where < zi < Z2 < ■ . . < Z2n- Denote by xi, . . . ,xn the odd labelled particles with 
respect to this ordering, and by yi,...,yiy the even labelled particles with respect to 
this ordering. We will also denote the probability density function for this ensemble by 
p^"'\xi,...,XN;yi,---,yN)- Thus 



2N 



p(")(a;i, . . . ,XAr;?/i, . . . ,?/Ar) = const- JJ 



e 2 z- 



i=l 



l<j<k<2N 



The {ki, /c2)-point correlation function for ki odd labelled particles, and k2 of even-labelled 
particles is defined as 



(iV-fci)! {N-k2)\ 



N N 



J(0,+oo)^-i J{0,oo)'=2 i=fcl+l.=fc2+l 



As it is explained in Forrester and Rains [22], p^^^ ^^^(xi, . . . , Xk^; . . . , 7/^2) describes 
parity respecting correlations of particles from the Laguerre orthogonal ensemble, see [22] 
for details and further references. Let P(kiM)i^^^ • • • ' ^^i' yii ■ ■ ■ 1 yk2) be similarly defined 
correlation function for the Meixner orthogonal ensemble. Let c — > 1^ and assume that 



Xi 




•^ki 


_l-c_ 


, . . . , 


_l-c_ 



are even, 



yi 




" yk2 


_l-c_ 


, . . . , 


l-c_ 



are odd. 
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Then definitions of tlie Meixner and Laguerrc orthogonal ensembles imply 



Theorem 17.1. We have 



1-c 



1-C 



yi 



l-c 



yk2 

1 - c 



PikiM)^^i^---^^k,;yi,---,yk2) 



Pf 



Kmiyj^^i) 



j,l=l,...,ki 
j=l,...,fc2;i=l,...,fci 



Kmiyj^yi) 



J i=l,...,A:i;Z=l,...,fc2 



j,l=l,...,k2 



where 



ss 



ATI 



K^Nli^,y) 



Ki^li^,y) 



ss 



(a) 
Nl 



SS 



(a) 
Nl 



SS'j^lS - s 



ss 



(a) 
Nl 



ss';^>s - s 



ix, y) 

ee 

ix,y) 

x, y) 

oe 

{x,y) 

x,y) 

eo 

{x,y) 
x,y) 

oo 

ix,y) 



c(") 

'-'ATI 

c(") C 
iJNI^ 



•^Nl 
c(") c 

cia) 
'-'Nl 

q(«) 

•^Nl 
c(") c 



{x,y) 

e 

ix,y) 
{x,y) 

e 

ix,y) 

{x,y) 

o 

{x,y) 
ix,y) 

o 

ix,y) 



and 



^^(a) 
SS 



Nl 

(a)" 
Nl 

C'(a) c 

iJni^ 

'JNI'^ 



SS 



(a) 



SS 



Nl 

(a) 
Nl 

qia) p 
iJNI^ 

q{a) p 

Jni'^ 



{x,y) = K^^\x,y) + ^/2N{2N + a) (f^^^^) {x)4''\y) = 
{x,y) = K^^\x,y) + ^2N{2N + a) = 
{x,y) = K^^\x,y) + ^2N{2N + a) = 

(x,y)-£:^(x,y)+V2iV(27V + a) (f^^'f ^) (x) (s^i^t^) (y), 
(a;,|/)-£^(a;,|/) + V27V(2Ar + a) (x) (y), 



(a;, 2/), 
{x,y), 



^i^Nl^ 


-s 


ee 

ix,y) = 


ce T^{a) 


c q{a) c 


-s 


eo 

{x,y) = 


cejy-ia) 


c c(a) c 
^^Nl^ 


-s 


oe 

{x,y) = 


cOT^ia) 
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(a) 
N 



s 



(a) 



S' 



Nl 
-1 ee 



(a) 
Nl 



{x,y) = 



x,y)-£''ix,y) + ^2Ni2N + a) (x) (y), 



5' 



iVl 



{x,y) 



q{a) 



[x,y) 



S 



(a) 
Nl 



ix,y). 



In the formulae above the operators £^ and £° are defined by the relations 



+ 00 



fiy)dy, {£''f){x) 



f{y)dy. 



Proof. The correlation kernel K^i for the Meixner orthogonal ensemble is obtained ex- 
plicitly in Theorem 13.81 a). To find the kernel for the correlation function ^^-^ we need 
to compute K^i{\yz^, \i^) as c — 1~. This can be done exploiting formulae in Lemma 
116.21 and in addition the following limiting relations 



lim —= 



lim 



(eK 



N 



X 



X 



1 - c 



{x) 



IS even, 



£y. 



2 



[X 



£''K 
£°K 



N 
N 



is odd, 

ix,y), [y^] is odd. 



IS even, 



□ 



Observe that the structure of the parity respecting correlation functions in Theorem 
117.11 is very similar to that in Section 3.1 of [22]. For a = the kernels of Theorem 117.11 
and the case A = of [22] must be equivalent. However, direct verification of this fact 
is not an easy task, see e.g. Section 5.1 of [22], and we postpone the discussion of this 
equivalence until a later publication. 
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